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Abstract. An extensive radiative transfer analysis of circumstellar SiO 'thermal' radio line emission from a large 
sample of M-type AGB stars has been performed. The sample contains 18 irregulars of type Lb (IRV), 7 and 
34 semiregulars of type SRa and SRb (SRV), respectively, and 12 Miras. New observational data, which contain 
spectra of several ground vibrational state SiO rotational lines, are presented. The detection rate was about 60% 
(44% for the IRVs, and 68% for the SRVs). SiO fractional abundances have been determined through radiative 
transfer modelling. The abundance distribution of the IRV/SRV sample has a median value of 6xl0~^, and a 
minimum of 2x10"^ and a maximum of 5xlO"^ The high mass-loss rate Miras have a much lower median 
abundance, <10~^. The derived SiO abundances are in all cases well below the abundance expected from stellar 
atmosphere equilibrium chemistry, on average by a factor of ten. In addition, there is a trend of decreasing 
SiO abundance with increasing mass-loss rate. This is interpreted in terms of depletion of SiO molecules by the 
formation of silicate grains in the circumstellar envelopes, with an efficiency which is high already at low mass-loss 
rates and which increases with the mass-loss rate. The high mass-loss rate Miras appear to have a bimodal SiO 
abundance distribution, a low abundance group (on average 4xl0~^) and a high abundance group (on average 
5xl0~^). The estimated SiO envelope sizes agree well with the estimated SiO photodissociation radii using an 
unshielded photodissociation rate of 2.5 x 10^^" s~^. The SiO and CO radio line profiles differ in shape. In general, 
the SiO line profiles are narrower than the CO line profiles, but they have low-intensity wings which cover the full 
velocity range of the CO line profile. This is interpreted as partly an effect of selfabsorption in the SiO lines, and 
partly (as has been done also by others) as due to the influence of gas acceleration in the region which produces 
a significant fraction of the SiO line emission. Finally, a number of sources which have peculiar CO line profiles 
are discussed from the point of view of their SiO line properties. 
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1. Introduction 

The atmospheres of and the circumstellar envelopes 
(CSEs) around Asymptotic Giant Branch (AGB) stars 
are regions where many different molecular species and 
dust grains form efficiently. The molecular and grain type 
setups are to a large extent determined by the C/0- 
ratio of the central star. For instance, SiO is formed 
in the extended atmospheres of both M-type [C/0<1; 
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* Based on observations using the SEST at La Silla, Chile, 
the 20 m telescope at Onsala Space Observatory, Sweden, the 
JCMT on Hawaii, and the IRAM 30 m telescope at Pico Veleta, 
Spain. 



0-rich] and C-type [C/0>1] AGB stars, but its abun- 
dance is much higher in the former. Therefore, the SiO 
'thermal' line emission (i.e., rotational lines in the v—0 
state; the term 'thermal' is used here to distinguish the 
v=0 state emission from the strong maser line emission 
from vibrationally excited states) is particularly strong 
towards M-stars, with the intensity of e.g. the J=2^1 
line comparable to, or even stronger than, that of the 
CO J=1^0 emission. Nevertheless, the initial observa- 
tions of SiO thermal radio line emission from AGB-CSEs 
llLambert & Vanden Bouil Il978t IWolff fc Carlsoni Il982h 
and their interpretation ijMorris et alJ . Il979j) suggested 
circumstellar SiO abundances (several) orders of magni- 
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tude lower th an tho se expected from the chemical equilib- 
rium models l|Tsuiil . [l973)) . 

Over the years the observational basis has 
proved consider ably 



, -Buiarrabal et al 

Bieging fc Latteil llQM iBieging etal 



im- 
1989^ 
1998L I200C 



1986, 



.01ofespn 'et'ldT n998ll. and even some inter f erome ter data 
exist ()Lucas et all Il993 ISahai fc BiegingL Il993|) . These 
data suggest that the SiO line emission originates in two 
regions, one close to the star with a high SiO abundance, 
and one extended region with a low SiO abundance. The 
relative contributions to the SiO line emission from these 
two regions depend on the mass-loss rate. 

This structure has been int erpreted as due to ac- 
cretion of SiO o n to du st grains l|Buiarrabal et all Il989l: 
ISahai fc BiegingL Il993|) . After the grains nucleate near 
the stars, they grow in part because of adsorption of 
gas-phase species. In 0-rich CSEs, refractory elements 
like Si, together with O, are very likely the main con- 
stituents of the grains, which are identified through the 9 
and 18 /xm silicate f e atures in the infrared spectr a of th e 
stars ((Forrest et al 1 Il975t IPegourie fc Panoula.rl |l985). 
Therefore, molecules like SiO are expected to be easily 
incorporated into the dust grains. As a result, the SiO gas 
phase abundance should fall off with increasing distance 
from the star as SiO molecules in the outflowing stellar 
wind are incorporated into the grains. The depletion pro- 
cess is, however, quite uncertain since it does not proceed 
at thermal equilibrium. Eventually, photodissociation de- 
stroys all of the remaining SiO molecules. 

The grain formation is important not only for the 
chemical composition of the CSE, but also because it af- 
fects its dynamical state [the radiation pressure acts on 
the gr a ins wh ich are dynamically coupled to the gas, e.g., 
iKwokl lll97,ql) ]. The SiO radio line profiles are narrower 
than those of CO and have mostly Gaussian- like shapes 
(e.g., iBuiarrabal et al.l l|l98(l Il989l) l a fact suggesting 
that the SiO line emission stems from the inner regions 
of the CSEs, where grain formation is not yet complete 
and where the stellar wind has not reached its terminal 
expansion velocity. This result is corroborated by interfer- 
ometric observations which show that the size of the SiO 
line er aitting region is independent of the line-of-sight ve- 
locity I Lucas et alin"9921 Lucas et al. explained this as a 
result of a rather exten ded acceleration region. However, 
ISahai fc BiegiMllll993 '). using a more detailed modelling, 
were able to explain both the line profiles and the bright- 
ness distributions with a 'normal' CSE, i.e., with a rather 
high initial acceleration. 

Therefore, 'thermal' SiO radio line emission is a use- 
ful probe of the formation and evolution of dust grains in 
CSEs, a complex phenomenon that is yet not fully under- 
stood, as well as the CSE dynamics. 

In this paper we present a detailed study of SiO ra- 
dio line emission from the CSEs of a sample of M-type 
AGB stars. The sample includes irregular (IRVs), semireg- 
ular (SRVs) and Mira (M) variables. The IRVs and SRVs 
have alread y been studied in c ircumstellar CO radio line 
emission (jOlofsson et allEoQ2i) . yielding estimates of the 



stellar mass-loss rates. Using these estimates a radiative 
transfer modelling of the SiO radio line emission is per- 
formed. A complete analysis of the circumstellar CO and 
SiO line emission is done for the Mira sub-sample. 

2. Observations of the IRV/SRV sample 

2.1. The IRV/SRV sample 

The sample contains all the M-type IRVs and SRVs 
det ected in circumstellar CO radio line emission 
by i Kerschbaum fc Olofsson ' (1999) and lOlofsson et ail 
(120021). The original source s e lectio n criteria are described 
in IKerschbaum fc OlofssonI ()l999j) . but basically these 
stars are the brightest 60 ^m-sources (IRAS Sqo typ- 
ically above 3Jy, with IRAS quality flag 3 in the 12, 
25, and 60 /im bands) that appear as IRVs or SRVs 
in the G e neral Catalogue of Variable Stars [GCVS4; 
iKholopovl (|l990|) ]. The detecti on rate of circumstellar CO 
was rather high, about 60% (|01ofsson et all (|2002fl : 69 
stars detected ). The basic pr o pert ies of the stars are listed 
in .Kerschba^im fc OlofssonI ((l999il and lOlofsson et alJ 
(|2002l) . 

The distances, presented in Tabled) were derived using 
an assumed bolometric luminosity of 4000 L© for all stars. 
We are aware of the fact that such a distance estimate have 
a rather large uncertainty for an individual object but it 
is adequate fo r a statistical study o f a sample of stars (see 
discussion bv lOlofsson et all ((20021) V The apparent bolo- 
metric fluxes were obtained by integrating the spectral 
energy distributions ranging from t he visual data over the 
near- infrared to the IRAS-range ((Kerschbaum fc HronL 
Il996l) . 

2.2. The observing runs 

The SiO {v=0, J=2^1; hereafter all SiO transitions are in 
the ground vibrational state) data were obtained using the 
20 m telescope at Onsala Space Observatory (OSO) and 
the 15 m Swedish-ESO Submillimetre Telescope (SEST) 
on La Silla, Chile. At SEST, a sizable fraction of the 
stars were observed also in the SiO J=3^2 line, and 
four additional sources were observed with the IRAM 30 m 
telescope at Pico Veleta, Spain, in this line. The higher- 
frequency lines, J=5— >4 and 6— >5, were observed towards 
10 and 3 stars, respectively. The observing runs at OSO 
were made over the years 1993 to 2000, at SEST over the 
years 1992 to 2003, and at IRAM between October 18 
and 22 in 1997. Telescope and receiver data are given in 
Table ^ Tree and ?7nib stand for the representative noise 
temperature of the receiver (SSB) and the main beam ef- 
ficiency of the telescope, respectively. 

Two fiherbanks at OSO (256 x250 kHz, and 
512 x1MHz), two acousto-optical spectrometers at 
SEST (86 MHz bandwidth with 43 kHz channel sep- 
aration, and 1 GHz bandwidth with 0.7 MHz channel 
separation) , and a 1 MHz filter bank at IRAM were used 
as spectrometers. Dual beam switching (beam throws of 
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Table 1. Data on telescopes and receivers 



Telescope 


Frequency 


Beamwidth 








[MHz] 


["] 


[K] 




OSO 


86847 


42 


150 


0.55 


SEST 


86847 


57 


100 


0.75 


SEST 


130269 


39 


120 


0.65 


IRAM 


130269 


18 


150 


0.58 


SEST 


217105 


25 


600 


0.55 


SEST 


260518 


21 


800 


0.45 



about 11'), in which the source was placed alternately in 
the two beams, was used to eliminate baseline ripples at 
OSO and SEST, while a wobbler switching with a throw 
of 150" in azimuth was used at IRAM. Pointing and fo- 
cussing were checked every few hours. The line intensities 
are given in the main beam brightness temperature scale 
(T,„b), i.e., the antenna temperature has been corrected 
for the atmospheric attenuation (using the chopper wheel 
method) and divided by the main beam efficiency. 

2.3. Observational results 

A total of 60 stars were observed in circumstellar SiO line 
emission (i.e., about 85% of the stars detected in circum- 
stellar CO): 34 stars were detected in the SiO J=2^1 line, 
21 in the ,7=3^2 line, and 3 in the J=5-^4 and J=6— >5 
lines. Clear detections of SiO lines were obtained towards 
36 sources, i.e., the detection rate was about 60%: 8 IRVs 
(detection rate 44%) and 28 S RVs (detection rate 68%) 
were detected. Tables lA.ll and IA.2I in the Appendix list 
all our SiO observations. The names in the GCVS4 and 
the IRAS-PSC are given. The first letter of the code de- 
notes the observatory (IRAM, OSO, or SEST), the rest 
the transition observed. Another code reflects the 'success' 
of the observation (Detection, Non-detection). 

The stellar velocity is given with respect to the helio- 
centric (uhoi) and LSR frame [ulsr! the Local Standard of 
Rest is defined using the standard solar motion (B1950.0): 
U0=2Okms-\ 00=270.5°, <50=-h3O°]. The stellar veloc- 
ity, the expansion velocity, and the main beam bright- 
ness temperature were obtained by fitting the function 
^mbll — {{v* — Vz)/vc)'^]'^ to the line profile. The inte- 
grated intensity, /= jT^ibdv, is obtained by integrating 
the line intensities over the line profile. The uncertainty 
in / varies with the S/N-ratio, but we estimate that it is 
on avarage <15%. To this should be added an estimated 
uncertainty in the absolute calibration of about 20%. For 
a non-detection an upper limit to / is estimated by mea- 
suring the peak-to-peak noise (Tpp) of the spectrum with 
a velocity resolution reduced to 15kms^^ and calculat- 
ing /=15Tpp. The Q-column gives a quality ranking: 5 
(not detected), 4 (detection with very low S/N-ratio <3), 
3 (detection, low S/N-ratio «5), 2 (detection, good S/N- 
ratio « 10), and 1 (detection, very good S/N-ratio >15). 
Finally, in cases of complex velocity profiles the measured 



component is indicated in the form b=broad, n=narrow, 
b-|-n=total. 

All the spectra are shown in Figs. IB. II to IB. 41 The 

velocity scale is given in the heliocentric system. The ve- 
locity resolution is reduced to 0.5kms~^, except for some 
low S/N-ratio spectra where a resolution of lkms~^, or 
even 2 kms""'^, is used, and for some low expansion velocity 
sources for which 0.25 km s~^ is used. 

3. The Mira sample 

In order to make a more extensive study of circumstellar 
SiO line emission in the CSEs of M-type AGB-stars, a 
sample of 12 Mira variables with higher mass-loss rates 
was added. The distan ces are obtained u sing t he period- 
luminosity relation of IWhiteloc.k et'ZI (jn^. Through 
modelling of their circumstellar CO radio line emission 
(Sect. EJ, we determined that 4 of the Miras have very 
high mass-loss rates (> 10"^ Mq yr~^), 6 are intermediate 
to high mass loss rate objects (> 10~^ M0 yr~^) and 2 are 
low mass-loss rate sources (a few lO^^M0yr^^). 

For this sample data has been gathered from a num- 
ber of sources . The C0( J=1^0) data were taken from 
lOlofsson et J] lll998l) . while the C0(J=2^1, J=3^2, 
and J=4^3) data were obtained from the archive of the 
James Clerk Maxwell Telescope on Mauna Kea, Hawaii. 
The JCMT data are taken at face value after convert- 
ing to the main beam brightness scale. However, in the 
cases where there are more than one observation available, 
the derived line intensities ar e generally consistent with in 
±20% (as was found also bv lSchoier fc Olofe son ('200l[)) 
The SiO( J=2-^l) data were obtained from lOlofsson et alJ 
(1998). The SiO J=5^4 Une was observed in four objects, 
and the J=6— >5 line in one object using SEST with the 
same observational equipment and procedure as described 
above. 

The relevant observational results are summarized in 
Table 1X31 and the SiO spectra are shown in Fig. IB. 51 The 
names in the GCVS4 and the IRAS-PSC are given. The 
first letter of the code denotes the observatory (JCMT, 
OSO, or SEST), the rest the transition observed. 

4. Modelling of circumstellar line emission 

Apart from presenting new observational results on ther- 
mal SiO radio line emission from AGB-CSEs a rather 
detailed modelling of the emission will be performed. In 
some senses this is a more difficult enterprise than the 
CO line modelling. The SiO line emission predominantly 
comes from a region closer to the star than does the CO 
line emission, and this is a region where the observational 
constraints are poor. The SiO excitation is also normally 
far from thermal equilibrium with the gas kinetic temper- 
ature, and radiative excitation plays a larger role (hence 
the term 'thermal' is really not appropriate) . Finally, there 
exists no detailed chemical model for calculating the ra- 
dial SiO abundance distribution. These effects make the 
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Fig. 1. The line intensity ratio 

/(SiO,J=2^1)//(CO,J=1^0) as a function of a 
mass-loss rate measure [detections are shown as filled 
symbols, while tentati ve and negative results are shown 
as open symbols; see lOlofsson et all lll998l) l for more 
details] 



SiO line modelling much more uncertain, and dependent 
on a number of assumptions. 

The aim is to investigate to what extent the ther- 
mal SiO line emission is a useful probe of e.g. the dust 
formation and the CSE dynamics. There are observa- 
tional indications that this is the case but the interpre- 
tation is normally n ot straightforward. As an example, 
lOlofsson et al.l l)l998l) found that the line intensity ratio 
/(SiO, J=2^1)//(C0, J=l— >0) decreases markedly as a 
function of a mass-loss rate measure. Their results are re- 
produced here, but now including all stars of our IRV/SRV 
and Mira samples, Fig.^ A straightforward interpretation 
would be that the SiO abundance decreases with mass-loss 
rate due to increased depletion efficiency and hence this 
limits severely the SiO line strength. However, excitation 
may play an important role here, both for SiO and CO, 
and a detailed modelling is required. 

4.1. The method 

In order to model the circumstellar SiO line emission a 
non-LTE radiative tran sfer code based on the Monte Carlo 
method has been used jBfirned . It has been previ- 

ously used to model circ umstellar CO radio l ine emission 
in samples of both C- llSchoier fc OlofssonL I2OO0I I2OOII: 
ISchoier et all l2002(l and 0-rich l)01ofsson et all l2002() 
AGB-CSEs, and also to model the HCN and CN hne 
emission from a limit ed number of C-rich AGB-CSEs 
l|Lindavist et all l200fl() . 

4.2. The SiO molecule 

In the excitation analysis of SiO 50 rotational levels 
in both the ground and the first excited vibrational 
state are considered. The energy levels of this linear ro- 
tor are calculated using the molecular constants from 



iMollaaehababa et alJ ijlQQlf) . The radia tive rates are cal- 
culate d using the dipole moment from iRavmonda et alJ 
Collisiona l deex citation rates have been calculated 
by Turner et al.l l)l992fl 'l in the temperature range 20- 
300 K and up to J=20. The original data set has been 
extrapolated in temperature and to include levels up to 
J— 50 (Schoier et al., in prep.). 

4.3. The circumstellar model 

The CSEs around AGB-stars are intricate systems where 
an interplay between different chemical and physical pro- 
cesses takes place. This makes the modelling of circum- 
stellar radio line emission a quite elaborate task. In the 
analysis presented here, a relatively simple, yet realistic, 
model for the geometry and kinematics of the CSEs has 
been adopted. Below follows a short description of the 
main featu res of the circumstellar mo del. For more details 
we ref er to ISchoier fc OlofssonI l)200l|) and lOlofsson et alJ 

iooi). 

A spherically symmetric geometry of the CSE is 
adopted. The mass loss is assumed to be isotropic and con- 
stant with time. The gas expansion velocity is assumed to 
be constant with radius. There is a possibility that neither 
the mass-loss rate nor the expansion velocity are constant 
in the regions of interest here. This should be kept in mind 
when interpreting the results. There is growing evidence 
for mass-loss m odulations of AGB-stars on a time scale of 
abou t 1000 vr JMauron fc HugginsL l200nHMarengo et all 
I2OOII iFong et al^l200.1^ ■ and the CO line emission comes 
from a much larger region than that of the SiO lines, and 
hence averages over a longer time span. Furthermore, the 
SiO line emission comes from the inner part of the CSE, 
where it is likely that the gas has not fully reached the 
terminal velocity. We have not allowed for the presence 
of gas acceleration nor a time-variable mass loss in the 
modelling in order to limit the number of free parameters. 

The inner boundary of the CSE was set to 1 x 10^^ cm 
(~3R*). This parameter is specially important in the case 
of SiO where radiative excitation is expected to play a role. 
A turbulent velocity of 0.5kms~^ is assumed throughout 
the entire CSE (see discussion by Olofsson et al, (200^ ). 
The outer boundaries of the molecular abundance distri- 
butions are, for both CO and SiO, determined by pho- 
todissociation due to the interstellar UV radiation field. 
For CO we use the modelling of lMamon et all l|l988|) . The 
procedure for SiO is presented in Sect.El 

The radiation field is provided by two sources. The cen- 
tral radiation emanates from the star. This radiation was 
estimated from a fit to the spectral energy distribution 
(SED) by assuming two blackbodies, one representing the 
direc t stellar radiation and one t he dust-processed radia- 
tion teersch baum fc HronL 11996V In the case of optically 
thin dust CSEs the stellar blackbody temperature derived 
in this manner is generally about 500 K lower than the ef- 
fective temperature of the star. The dust mass-loss rates 
of the IRV/SRVs are low enough that the dust blackbody 
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can be ignored. For the sample of Mira variables both 
blackbodies were used, since for these high mass-loss rate 
stars, the excitation of the SiO molecules may be affected 
by dust emission. The second radiation field is provided 
by the cosmic microwave bakground radiation at 2.7 K. 

In the SiO line modelling the gas kinetic temperature 
law derived in the modelling of the circumstellar CO radio 
line emission was used. This is reasonable since the SiO 
line emission contributes very little to the cooling of the 
gas. However, the SiO line emission comes mainly from 
the inner CSE, where the CO lines do not put strong con- 
straints on the temperature, and where other coolants, 
specifically H2O, may be important. We estimate though 
that the kinetic tempartures used in our modelling are 
not seriously wrong. In addition, for at least the lower 
mass-loss rates the SiO molecule is mainly radiatively ex- 
cited, and hence the exact gas kinetic temperature law 
and the coUisional rate coefficients play only a minor role, 
see Sect. 14.41 

In Sect. 14.41 some implications of these assumptions are 
discussed. 

The mass- loss rates for the sample of IRV/SRVs were 
already presented in Olofsson et al. ( 2002). They were de- 
rived through modelling of circumstellar CO radio line 
observations. A median mass loss rate of 2x 10~^ Mq yr"-'^ 
was found for this sample. These mass-loss rate estimates 
are expected to be accurate to within a factor of a few 
for an individual object. Nevertheless, they are probably 
the best mass-loss rate estimates for these types of ob- 
iec ts [also in agreemen t with the mass-loss rate estimates 
by iKnapp et all lll99j^ for five sources m common . 

The modelling of the circumstellar CO radio line emis- 
sion for the Mira sa mple is presented in this paper. The 
same approach as in lOlofsson et al.l 1 20021) has been used, 
i.e., the energy balance equation is solved simultaneously 
with the CO excitation. A number of (uncertain) parame- 
ters describing the dust are introduced. They are grouped 
in a global parameter, the /i-parameter, which is given by 



0.01 



2.0 gem 

Pgr 



0.05 yum 



(1) 



where ip is the dust-to-gas mass ratio, pgr the dust grain 
density, and Ogr its radius. This parameter is particu- 
larly important for the heating due to gas-grain colli- 
sions. The normalized values are the ones used to fit 
the C O radio line emission of IRC+10216 using this 
model l^^hi^i^r ^ Olofsson* (20Q1D, i.e., h^l for this ob- 
ject. ISchoier fc Olof sson ( 200^ found that on average 
h=0.2 for the lower luminosity sources (below 6000 Lq; 
and h=0.5 for the more luminous sources) i n their sam- 
ple of bright carbon stars and Olofsson et alJ l|2002') found 
h=0.2 for their sample o f M-typ e IRV/SRVs. In addition, 
following lOlofsson et alJ l)2002j) we use (in the gas-grain 
drift heating term) a flux-averaged momentum transfer 
efficiency from dust to gas, Qp.F, equal to 0.03 indepen- 
dent of the mass-loss rate, and adopt a CO abundance 
with respect to H2 of 2 x 10"''. The latter may very well 



be an underestimate for these high mass-loss rate objects 
(see below). 



4.4. Dependence on parameters 

A sensitivity test has been performed in order to deter- 
mine the dependence of the calculated SiO line intensities 
on the assumed parameters for a set of model stars. They 
are chosen such that they have nominal mass-loss rate and 
gas expansion velocity combinations which are character- 
istic of our samples: a low mass-loss rate (10^^ M0 yr^^, 
7knis^^), an intermediate mass-loss rate (10~^MQyr~^, 
lOkms^^), and a high mass-loss rate (10~^MQyr~^, 
15kms~"'^) model star. They are placed at a distance of 
250 pc (a typical distance of the stars in the IRV/SRV 
sample). We have also taken nominal values for the lumi- 
nosity (i=4000 L0 for the low and intermediate mass- loss 
rate model stars, and i=8OOOL0 for the high mass-loss 
rate model star), the effective temperature (rbb=2500K), 
the /i-parameter {h=0.2 for the low mass-loss rate model 
star, and h—0.5 for the other two), the envelope inner ra- 
dius {ri=2 X lO^** cm, which is twice the inner radius used 
in the modelling), the turbulent velocity (ui=0.5kms~^), 
and the SiO abundance [/sio=5 x 10"^ (close to the me- 
dian value for our IRV/SRV sample, see below); through- 
out this paper the term abundance means the fractional 
abundance with respect to II2, the dominating molecu- 
lar species in the CSEs]. The SiO envelope outer radius 
is calculated for each model star following the same rela- 
tion that is used in the modelling of the sample stars (see 
Sect. 16.41 . The SiO lines are observed with beam widths 
characteristic of our observations. All parameters (except 
the mass-loss rate and expansion velocity) are changed by 
—50% and -1-100% and the velocity-integrated line inten- 
sities are calculated. In order to check the effect of the 
/i-parameter on the modelled intensities the radial gas ki- 
netic temperature law is scaled by —33% and -f 50%. The 
results are summarized in Table |21 in terms of percentage 
changes. To see how the SiO/CO line intensity ratios vary 
with mass-loss rate, the CO line intensities derived from 
the models with the nominal parameters are also included. 

Despite the fact that the dependences are somewhat 
complicated there are some general trends. The line in- 
tensities are, in general, sensitive to changes in the outer 
radius, but less so for the high-J lines, a fact which is 
more evident for the low mass-loss rate stars. There is 
also a dependence of all line intensities on the SiO abun- 
dance, irrespective of the magnitude of the mass-loss rate. 
These particular dependences of the line intensities on the 
envelope outer radius and the SiO abundance allowed us 
to derive envelope sizes for those stars with multi-line ob- 
servations (see Sect. 16.31 . The line intensities are rather 
insensitive to a change in the kinetic temperature. Only 
the high-J lines for high mass-loss rates show a weak de- 
pendence on this parameter. The dependence on the inner 
radius is marginal, and so is the dependence on the tur- 
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Table 2. The effect on the integrated model SiO intensities (in percent), due to changes in various parameters. Three 
model stars with mass loss rate and gas expansion velocity characteristics typical for our samples are used. They lie at a 
distance of 250 pc, and have luminosities of 4000 Lq (the model stars with mass-loss rates of 10"'' and 10~^ Mq yr~^) 
and 8000 Lq (the model star with a mass-loss rate of 10~^ Mq yr~^), and blackbody temperatures of 2500K. The 
nominal CSE parameters are h—0.2 (for the lowest mass-loss rate model star) and h—0.5 (for the other two model 
stars), ^=2x10^"* cm, z;t=0.5kms~-^, and /sio=5xlO~^. The size of the SiO envelope, r^, is given by Ea. ((T^ for the 
given mass-loss rate. The SiO J=2^1, J=3^2, J=5— >4, and J=6— >5 lines are observed with beam widths of 57", 38", 
23", and 19", respectively (appropriate for a 15m telescope). The model integrated line intensities, I in Kkms~^, are 
given for the nominal parameters. For comparison, also the integrated CO line intensities for the model stars are given 
[for a 15m telescope; 45" (J=1^0; 7 = 0.14, 3.8, and 45Kkms-i for lO-^lO'^, and lO-^Moyr"!, respectively), 23" 



(J=2- 


^1), 15" (J: 


=3^2), 9" (J 


=5- 


4), 7" 5 (J= 


=6^ 


5)] 














Par. 




10-^ 


Mo yr-^ 


, 7kms 


10-^ 


' Mo yr- 


\ lOkms"^ 


10"^ 


Mo yr 


\ 15kms"^ 


Change 


2-1 


3-2 


5-4 


6-5 


2-1 


3-2 


5-4 


6-5 


2-1 


3-2 


5-4 


6-5 


IsiO 




0.11 


0.30 


0.71 


0.90 


1.3 


2.7 


5.5 


6.8 


10 


20 


39 


48 


I CO 




1.6 


3.9 


6.0 


6.4 


17 


31 


46 


50 


120 


184 


252 


271 


fsiO 


-50% 


-40 


-37 


-36 


-36 


-33 


-30 


-31 


-32 


-27 


-24 


-26 


-27 




+ 100% 


+70 


+56 


+44 


+44 


+46 


+39 


+40 


+42 


+29 


+27 


+36 


+37 


L 


-50% 





-7 


-14 


-14 


-8 


-5 


-4 


-4 


-10 


-2 










+100% 





+11 


+ 19 


+21 


+13 


+14 


+ 10 


+9 


+19 


+8 


+4 


+3 




-33% 








-7 


-8 


-5 


-9 


-16 


-19 


-16 


-21 


-27 


-31 




+50% 





+4 


-12 


+7 


+8 


+11 


+16 


+18 


+20 


+24 


+32 


+35 


Tc 


-50% 


-60 


-56 


-39 


-33 


-48 


-38 


-29 


-25 


-35 


-33 


-26 


-22 




-hl00% 


+110 


+70 


+32 


+24 


+53 


+41 


+22 


+ 14 


+41 


+39 


+18 


+9 


fi 


-50% 


+10 


+4 


+2 


+3 


+2 


+2 


+ 1 





+ 1 


-1 


-1 







+100% 





-4 


-8 


-11 


-1 





-3 


-4 


+ 1 


+1 


-1 


-3 


Vt 


-50% 








-2 


-1 


+1 


+1 


-1 


+2 


-6 


-6 


-2 


-2 




+100% 





+4 


+2 


+3 





+3 


+3 


+2 


+3 


+3 


+ 1 






bulent velocity width (as long as it is significantly smaller 
than the expansion velocity). 

The dependece on luminosity is also weak, with only 
small changes in high- J line intensities for low mass-loss 
rates and in low- J line intensities for high mass-loss rates. 
However, the radiation field distribution may be of im- 
portance here, in particular for the high mass-loss rate 
objects. We have checked this for the high mass-loss rate 
model star. If half of the luminosity is put in a 750 K 
blackbody, the J=2^1, J=3^2, J=5^4, and J=6^5 
line intensities increase by a factor of 1.7, 1.3, 1.1, and 
1.1, respectively. That is, the lower J-lines are most af- 
fected, partly because of maser action (in particular in 
the J=1^0 line). This means that the SiO abundance 
estimates for the high mass-loss rate Miras are particu- 
larly uncertain, and the line saturation makes things even 
worse. 

A velocity gradient may affect the SiO line inten- 
sities since it allows the central pump photons to mi- 
grate further out in the CSE. We have tested a velocity 
law of the form (ap propriate for a dust-driven wind, see 
iHabing et all 

v{r) = ^vf + {v^^-vf)(^l^!l^ , (2) 

where Vi is the velocity at the inner radius, and v^c the 
terminal velocity. This produces a rather smooth increase 



in velocity, and for Ui/woo=0.25 (which we have used) 
90% of the terminal velocity is reached at r = lOn (for 
low mass-loss rate objects, this is also the region which 
produces the main part of the SiO radio line emission). 
There is only an effect for the low mass-loss rate object and 
the higher- J lines. For instance, the J— 6—^5 line intensity 
increases by about 10%. A velocity gradient has though 
the effect that the lines become narrower (Sect, mj) . 

Finally, the line intensity ratio /(SiO, J=2^1)/ 
/(CO,J=1-^0) decreases with mass-loss rate: 0.79 for a 
mass-loss rate of 10~'^MQyr-\ 0.33 for 10"^MQyr~\ 
and 0.24 for lO^'^Mgyr"!. This result is in line with the 
observational result presented in Fig.Q] and suggests that 
at least part of the trend is an excitation effect. 

5. CO modelling of the Miras 

In order to obtain mass-loss rates for the Mira sample we 
have modelled the circumstellar CO radio line emission ob- 
served toward s these stars u sing the procedure described 
above and in ISchdier fc Olofsson (.2001 ) . The estimated 
mass-loss rates are given in Table 13 rounded off to the 
number nearest to 1.0, 1.3, 1.5, 2.0, 2.5, 3, 4, 5, 6, or 8, 
i.e., these values are separated by about 25%. The dis- 
tribution of derived mass-loss rates have a median value 
of 1.3x 10~^ Mq yr~^. Therefore, these Miras sample the 
high mass-loss rate end of AGB stars. Only two of them 
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Source 


P 






T, 


Td 


D 




M 


'p 


h 


Xrcd 


N 




[days] 


[Lo] 




[K] 


[K] 


[pc] 


[km s-i] 


[10-s Mq yr-i] 


[10^*^ cm] 








TX Cam 


557 


8400 


0.26 


1800 


460 


380 


18.5 


7 


8.7 


1.5 


0.5 


4 


R Gas 


431 


6500 


0.06 


2100 


530 


220 


10.5 


1.3 


4.0 


1.1 


1.0 


4 


R Hya 


388 


5800 


0.03 


2300 


580 


150 


7.0 


0.3 


2.0 


0.5 




1 


R Leo 


313 


4600 


0.02 


2100 


570 


130 


6.0 


0.2 


1.7 


0.6 


0.4 


3 


GX Mon 


527 


8000 


0.38 


1500 


380 


540 


18.7 


40 


24.0 


0.5 


2.4 


4 


WX Psc 


660 


10000 


0.85 


920 


400 


600 


19.3 


10 


11.4 


0.4 


6.6 


4 


IK Tau 


500 


7500 


0.46 


1500 


500 


250 


18.5 


30 


20.5 


0.3 


0.4 


4 


IRG+ 10365 


500 


7500 


0.26 


1600 


430 


750 


16.2 


30 


23.8 


0.5 


2.2 


2 


IRC -10529 


680 


10400 


0.17 


1000 


410 


270 


12.0 


2.5 


5.8 


0.1 


1.6 


4 


IRC-30398 


575 


8700 


0.24 


2000 


480 


390 


16.0 


6 


8.2 


0.5 


0.1 


2 


IRC+40004 


750 


11500 


0.44 


1700 


400 


410 


18.0 


6 


8.6 


0.5 


2.4 


2 


IRC+50137 


629 


9600 


0.57 


1300 


310 


410 


17.0 


10 


10.7 


0.1 


4.2 


3 



^ Derived from a period-luminosity relation 
^ The GO photodissociation radius 



(R Hya and R Leo) have low to intermediate mass-loss 
rates (a few times lO~^M0yr~^). h was used as a free 
parameter in the fit for those sources with more than two 
lines obser ved. The average value is 0.6, i.e., very simi- 
lar to what lSchoier fc OlofssonI ll200lh found for the more 
luminous stars in the their carbon star sample. We used 
h=0.5 for those stars observed in only one or two lines. 
The quality of the fits are given by the chi-square statistic 
Xred (^^'^ Sect. l7.1l for the definition). 

A CO fractional a bundance of 2x l 0~"^ h as been used 
following the work of lOlofsson et alJ l)2002j) on the CO 
modelling of low to intermediate mass-loss rate IRV/SRVs 
of M-type. It is quite possible that, for the high mass-loss 
rate stars involved here, the CO abundance is higher due 
to a more efficient formation of CO at higher densities and 
lower temperatures. A higher CO abundance would lower 
somewhat the derived mass-loss rates. 

Among the Miras with the highest mass-loss rates 
there is a trend that the model J=1^0 line intensities 
are low for a model which fits well the higher- J lines. The 
reason is that the CO lines reach the saturation regime at 
about 10^^ Mq yr~^, with the higher-J lines saturating 
first. Therefore, we chose to put more weight on the high- 
J lines in the model fit. The reported values for the mass- 
loss rates of these stars are, in this context, therefore con- 
sidered to be lower li mits. This typ e of pr oblem has also 
been encountered by iKemner et all (j2Qo5). For WX Psc, 
the only star in common with us, they derived a mass-loss 
rate of 1.1x10^''' Mq yr^^ by fitting the J=2^1 line, and 
successively lower mass- loss rates for the higher-J lines 
reaching about 10~^ M0 yr"-'^ by fitting the J~6-^5 and 
J— 7^6 lines. In this work a value of 1.1 x 10^^ M0 yr^^ is 
derived based on the J=2— ^1, J— 3-^2, and J— 4^3 lines, 
but a fit to the J— 1^0 line requires a mass- loss rate about 
a factor of three higher. Kemper et al. speculate that vari- 
able mass loss and gradients in physical parameters (e.g., 
the turbulent velocity width) may play a role. To this we 



add that the size of the CO envelope, which mainly affects 
low-J lines, is important. 

The CO expansion velocities given in Table El are ob- 
tained in the model fits. Hence, they are somewhat more 
accurate than a pure line profile fit, since for instance the 
effect of turbulent broadening is taken into account. The 
uncertainty is estimated to be of the order ±10%. The 
gas expansion velocities have a distribution with a me- 
dian value of 15.3kms~^, while the IRV/SRV sample has 
a median gas expansion velocity of 7.0 km s~^ . Again, only 
R Hya and R Leo have low CO expansion velocities, below 
lOkms-^ 



6. Size of the SiO envelope 

The results of the SiO line modelling will depend strongly 
on the adopted sizes of the SiO envelopes. Unfortunately, 
these are not easily observationally determined nor the- 
oretically estimated. Early work assumed that the whole 
CSE contributes to the ob served SiO thermal line emission 
fe.g.- lMorris et al I ("1979^). The mostly Gauss an-like SiO 
profiles found by [Buiarrabal ct al. (1986, 1989) towards 



0-rich CSEs suggested that this is not the case. The gen- 
erally small size of the SiO thermal line emitting region 
requires interferometric observations in order to resolve it. 
Results from SiO multi-line modelling and interferometric 
data will be combined here to estimate the sizes of the 
SiO envelopes. 

6.1. The SiO abundance distribution 

Previous work strongly suggests that the SiO abundance 
in the CSE is markedly lower than that in the stellar at- 
mosphere. The decrease in the SiO abundance with ra- 
dius is very likely linked to two different processes taking 
place in the CSE. Photodissociation due to interstellar 
UV radiation is a well-known mechanism which reduces 
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the abundances of molecules in the extended CSE, but 
for SiO the depletion onto grains closer to the star must 
also be taken into account. We outline here in a simpli- 
fied wa^;_th£_effects of the se processes [based on th e works 
by I.Tura fc Morris ^98^; iHuggins fc Glassgoldl (|l982)]. 
However, the theoretical results are not used in our mod- 
elling, but they serve as a guide for the assumptions and 
the interpretation. 

Since the rate of evaporation is very large for 
Tgr > (Tbind/50) (where Tgr is the grain temperature, and 
fcTbind the binding energy of the molecule onto grains), 
there is a critical radius, rg, such that for smaller radii 
there is effectively no condensation, while for larger radii 
almost every molecule that sticks onto the grain remains 
there. The value of tq can be estimated from the condition 
that the characteristic flow time, r/uo(r), is equal to the 
evaporation time [i?evap(Tgr)]~^- A classical evap oration 
theor y has been used to obtain the rate for CO ijLegerl 
Il98,ll) . and the result is 



f = 



3xl0i3exp(-^ 



cm. 



(3) 



where Vc{ro) is given in cms^^. While different species 
have different coefficients in front of the exponential, by 
far the most important term is the exponential. The rate 
of classical evaporation is generally so large that unless 
TgT < 7bind/50, condensation onto grains is not impor- 
tant. Therefore, in describing the condensation process, 
only variations in Tbind for different substances are con- 
sidered and variations, among species, of the constant co- 
efficient in Eg . O l are ignored. Tbind = 29500 K for SiO 
i Leger et al.l Tg,{r) = T^{R^/2rf-* (appropriate 

for an optically thin dust CSE), and Wc^lOkms^^ re- 
sults in a typical condensation radius for our sources (with 
L = 4000Lq and rbb=2500K) of about SxlQi^cm. 

Using the formulation by Jura & Morris (1985), the 
radial variation of the SiO abundance in a CSE, taking 
into consideration the depletion of molecules onto dust 
grains, is given by 



/sio('') = /sio(''o) exp 



where Tgcaic is a scale length defined by 



' scale 



(4) 



(5) 



where a is the sticking probability of SiO onto grains, 
the dust mass-loss rate in terms of dust grain number, dgi. 
the grain cross section, and Vdr the drift velocity of the 
dust with respect to the gas, obtained from the formula 



Vdr 



LvcQ 



cM 



(6) 



Thus, the abundance decreases due to condensation until 
it reaches the terminal value 



/sio(oo) = /sio(-R*)e 



-facalo/l-O 



(7) 



The condensation efficiency depends strongly on the 
dust mass-loss rate. For instance, ip — 0.002 (appro- 
priate for the average ft.-value of the IRV/SRV sam- 
ple), agr = 0.05/im, pgr = 2 gcm^'^, Qp_F = 0.03, a—1, 
L = 4OOOL0, Tbb=2500K, and We^lOkms^^ result 
in /sio(oo)//sio(-R*)-values of 0.76, 0.42, and 0.07 
for mass loss rates of 10"'' yr~^, 10~^ Mq yr~^, 
and lO^^MQyr"^, respectively. The corresponding 
/sio(oo)//sio(^*)-values for ^ = 0.01 is 0.38, 0.013, and 
10~^. Thus, we expect condensation to play only a minor 
role for the low mass-loss rate objects, but its importance 
increases drastically with the mass-loss rate. 

The particular radius at which the photodisso- 
ciation becomes effective depends essentially on the 
amount of dust in the envelope, which provides shield- 
ing against the UV radiation, and the abundance of vari- 
ous molecular species if the dissociation occurs in lines. 
iHuggins fc Glassgoldl l)l982() describe the radial depen- 
dence of the abundance of a species of photospheric origin 
that is shielded by dust (in the case of SiO, shielding due to 
H2O may be important but we ignore this here). Adopting 
this description in the case of SiO the result is 



dfsio 
dr 



G, 



0,SiO 



exp 



/siO, 



(8) 



where fsio is the fractional abundance of SiO with re- 
spect to H2, Go, SiO the unshielded photodissociation rate 
of SiO, and da,,q the dust shielding distance given by (see 
lJura fc Mo^ ljl98l|) ) 



^SiO 



1.4 



3(Q/agr)siO Md hM 



4pg 



Attva 



oc 



^^d 



(9) 



where Q is the dust absorption efficiency, Afd the dust 
mass-loss rate, and Ud the dust expansion velocity given by 
Vc + Vdr- The abundance decreases roughly exponentially 
with radius and we adopt /sio(^p) — /sio(^^*)/c to define 
the outer radius rp. It is obtained by solving the equation 



0,SiO 



E2{dsio/rp) ' 



(10) 



where i?2(x) is the exponential integral. 

Most likely the radial distribution of the SiO molecules 
is determined by a combination of the condensation and 
photodissociation processes. Thus, one can imagine an 
initial SiO abundance determined by the stellar atmo- 
sphere chemistry. For low mass-loss rates, the abundance 
decreases only slowly beyond the condensation radius un- 
til the photodissociation effectively destroys all remaining 
SiO molecules. For high mass-loss rates, the abundance de- 
clines exponentially beyond the condensation radius with 
an e- folding radius that can be estimated from Eq. 



1 



1 



^0 ^scale 



(11) 



(applicable only when rgcaic > ''o)- Using the same param- 
eters as above, except ■0 = 0.005 (appropriate for the aver- 
age /i- value of the high mass- loss rate stars), L = 8000 Lq, 
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and Wo = 15kms^^, the result for lO^^M0yr~^ is 
rc = 10^^ cm, i.e., only about twice the condensation ra- 
dius. An abundance decrease by a factor of a hundred is 
reached at about 4x 10^^ cm, which is about a factor of five 
smaller than the estimated SiO photodissociation radius 
for such an object. Once again, the results are sensitively 
dependent on the dust parameters, e.g., V' = 0.002 results 
in = 2x 10^'^ cm, but the abundance (before photodisso- 
ciation) never decreases by more than a factor of five. 

6.2. The adopted SiO abundance distribution 

For the radial distribution of the SiO abundance in the 
CSEs we adopt a Gaussian fall-off with increasing distance 
from the star, 

/siO = /ce-('-/'^")', (12) 

where fc is the central abundance, and To the e-folding 
distance. 

This is a considerable simplification to the complicated 
SiO abundance distribution. However, as shown above, the 
expected distribution depends so sensitively on the pa- 
rameters adopted (in particular the dust mass loss rate) 
that a more sophisticated approach is, for the moment, 
not warranted. We expect Eq. I|12|l to be a reasonable 
approximation to the SiO abundance distribution inside 
the photodissociation radius for the low and intermediate 
mass-loss rate objects. Equation (|12() is a reasonable ap- 
proximation for also the high mass-loss rate objects, but 
the size is either determined by condensation (high ■0) or 
photodissociation (low -0). 

We have checked whether the region within the con- 
densation radius, with a high SiO abundance, contributes 
substantially to the observed line intensities. For the 
model stars used in Sect. 14.41 it is found that a high SiO 
abundance (5 x 10^^) inside the condensation radius con- 
tributes by at most 20% of the line intensities from the 
rest of the SiO envelope. 




[cm] 



Fig. 2. contours (at the 1, 2, and 3cr levels) in the 
SiO abundance and envelope size plane for two stars. In 
the case of RX Boo there is a sufficient number of lines, 
including high-J ones, to constrain the size of the SiO 
envelope. In the case of R Dor the high-J line data are 
missing and only a lower limit can be obtained (since we 
do not expect the SiO envelope to be larger than the CO 
envelope) 



6.3. Results from SiO line modelling 

The model code used in this work allow us to estimate 
SiO envelope sizes provided that multi-line SiO observa- 
tions are available. The emission from higher- J lines comes 
very likely from the warmer inner regions of the SiO enve- 
lope. Therefore, the intensities of these lines can be fitted 
by varying only the SiO abundance, i.e., they are rather 
insensitive to the outer radius of the SiO envelope (see 
Table El). Once the SiO abundance has been found, the 
lower- J lines can be used as constraints to derive the size 
of the SiO envelopes, since their emission is photodissoci- 
ation limited (i.e., not excitation limited). 

It turns out that high-J line data, e.g., J=8— >7, 
are required to constrain both the abundance and the 
size. These c r ucial high-J line data were taken from 
iBieging et all (|200nl) . In the case of data including only 
moderately high-J lines, e.g., J=5— >4, only a lower hmit 



to the size can be obtained. This is illustrated in Fig. [5] 
where maps are given for two cases (see the definition 
of the chi-square statistic below). In this way, through 
the use of maps, we managed to estimate the SiO 
envelope sizes in 4 cases (RX Boo, R Gas, IRG— 10529, 
IRG-t-50137), and obtain lower limits to them in 7 cases 
(TX Gam, R Grt, R Dor, R Leo, GX Mon, Pup, 
IRG-30398). 

The resulting r^'S from the modelling are plotted as a 
function of the density measure M/vc, in Fig.^ We have 
here chosen to use the lower limits to the SiO envelope 
sizes for all sources in order to be consistent. The mini- 
mum least-square correlation between these SiO envelope 
radii and the density measure is 

log To = 19.2 + 0.48 log ( — I , (13) 
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Fig. 3. SiO photodissociation radii (obtained using the 
unshielded photodissociation rate Go = 2.5 x 10^^" s~"^) 
versus the estimated sizes of the SiO envelopes (IRV: 
square; SRV: circle; Mira: triangle). The dashed line shows 
the one-to-one relation 



(the correlation coefficient is 0.83) where is given in cm, 
M in Mq yr"-'^, and in kms~^. For the rest of the sources 
the SiO envelope sizes could not be derived through mod- 
eUing. 

We have checked our model results against those of the 
photodissocation model. The photodisso ciation rad ii are 
estimated from Eq. Hl()|l assuming and 
using the appropriate M- and /i- values for each source. 
A very good agreement with the estimated SiO enve- 
lope sizes (for all sources with detected SiO lines), from 
Eq. (|13|l . is obtained with an unshielded photodissocia- 
tion rate Go,SiO = 2.5 x 10"-^° s^^ (the average deviation 
is about 30%), see Fig. El This value is low er by about a 
factor of two to three than those reported bv lvan Dishoeckl 
l)l988|) and [Tarafdar & Dalgarno (1993), and higher by 
about a factor of two than that reported bv lLe Te^iff et al] 
l|2000l) . The latter report an uncertainty by (at least) a 
factor of two in their estimate. Thus, within the consider- 
able uncertainties, our line modelling results are in excel- 
lent agreement with those of the photodissociation model. 
The Tpis for our sample are given in Tabled On average, 
the photodissociation radii of SiO are about a factor of 
6 smaller than those of CO (the CO results are given in 
citeolofetal02. 



6.4. Interferometry data 

We have also checked our modelling results by compar- 
ing with the interferomctric SiO(J=2— > 1) dat a toward a 
number of 0-r ich CSEs of Lucas et alJ lll992l) . They de- 
rived the sizes of the SiO line emitting region from direct 
fits, assuming exponential source-brightness distributions, 
to the visibility data. Their observations thus yielded the 



XI 




10 " 10 ' 10" 

M/Vg [Mq yr"^ km"^ s] 

Fig. 4. The sizes of the SiO envelopes estimated from the 
SiO line modelling are plotted versus a density measure 
(open circles). The dashed line gives the relation between 
the SiO envelope size and the density measure given in 
Eq. I|13|l . Half intensity radii derived from interferomctric 
SiO(J=2— >1) observations are shown as solid circles. The 
solid line is the fit to the data given in Eq. (|14|l 

half-int ensity angular radii of the SiO(J=2^1) emitting 
regions. ISahai fc Biegind l|l993l) observed a smaller sample 
of CSEs interferometrically, and claimed that the source 
brightness distribution is rather of a power-law form (i.e., 
scale- free). This would explain why Lucas et al. derived 
essentially the same angular sizes for most of the sources 
independent of their distances. To resolve this issue re- 
quires more detailed observations, and we will only use 
the results of Lucas et al. to compare with our modelling 
results. 

We have six stars in common with iLucas et al ] ill 9921) 

(RX Boo, R Cas, W Hya, R Leo, WX Psc, IK Tau). FigH 
shows the intensity radii as a function of the density mea- 
sure M/ve, using our derived mass-loss rates, gas expan- 
sion velocities, and distances. The minimum least-square 
correlation between these intensity radii and the density 



measure is 



logr/ya = 18.8 + 0.49 log [ — 



(14) 



(the correlation coefficient is 0.88) where r//2 is given in 
cm, M in MQyr~^, and Vc in kms""'^. 

Thus, the scaling with the density measure of the in- 
tensity radii is in perfect agreement with our modelling 
result for the envelope sizes. The estimated SiO envelope 
sizes that are required to model the data are about three 
times larger than the SiO( J=2^1) brightness region. This 
may at first seem somewhat surprising, but a test us- 
ing the lO~^M0yr~^ model star of Sect. 14.41 which has 
an SiO envelope radius of I'.'T, shows that the resulting 
SiO(J=2^1) brightness distribution has a half-intensity 
radius of 0'.'4, i.e., about four times smaller. 
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7. Results of the SiO line model fits 

7.1. The fitting procedure 

The radiative transfer analysis produces model brightness 
distributions. These are convolved with the appropriate 
beams to allow a direct comparison with the observed 
velocity-integrated line intensities and to search for the 
best fit model. As observational constraints we have used 
the data present ed in this paper and the high-frequency 
data obtained bv lBieeing et al.1 1)200(]|) . With the assump- 
tions made in the standard circumstellar model and the 
mass-loss rate and dust properties derived from the mod- 
elling of circumstellar CO emission, there remains only one 
free parameter, the SiO abundance [for all stars is taken 
from Eq. ifT^ ]. The SiO abundance was allowed to vary in 
steps of RilO% until the best-fit model was found. The 
quality of a particular model with respect to the observa- 
tional constraints can be quantified using the chi-square 
statistic, 



Xrcd 



1 



P 



N 

E 

1=1 



[In 



(15) 



where / is the total integrated line intensity, Ui the uncer- 
tainty in observation i, p the number of free parameters 
(2 in the cases of multi-line CO modelling, but only 1 for 
the SiO line modelling, except in the cases discussed above 
where also was a free parameter), and the summation 
is done over all independent observations N. The errors 
in the observed intensities are always larger than the cal- 
ibration uncertainty of «20%. We have chosen to adopt 
Ci = 0.2/obs,i to put equal weight on all lines, irrespective 
of the S/N-ratio. The final chi-square values for stars ob- 
served in more than one transition are given in Table ^ 
They are, in many cases, rather large suggesting that our 
circumstellar model may not be entirely appropriate for 
the modelling of the SiO radio line emission, see Sect. |H1 
The line profiles were not used to discriminate between 
models, but differences between model and observed line 
profiles are discussed in Sects 17.41 and ISl 

7.2. The accuracy of the estimated abundances 

We will here try to estimate the uncertainty in the derived 
SiO abundances. The uncertainties due to the adopted cir- 
cumstellar model are ignored since these are very difficult 
to estimate, and focus is put on those introduced by the 
adopted parameters (see Sect. 14.41 . We start by consid- 
ering the IRV/SRVs. The results depend crucially on the 
validity of Eq. lO. A change by -50% and -1-100% in 
the size of the SiO envelope results in a variation of the 
J=2^1 line intensity by about ±50%, and therefore an 
equal uncertainty in the abundance. The product of /gio 
and M is essentially constant for a best fit model. It is 
estimated that the mass-loss rate is uncertain by at least 
a factor of two (due to the modelling). An uncertainty in 
the distance has only a minor effect on the abundance (the 
change in mass-loss rate compensates for the change in 
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Fig. 5. SiO fractional abundances versus the mass- loss 
rate (IRV:square, SRV:circle, Miraitriangle). The horizon- 
tal line marks the maximum abundance allowed by solar 
abundances. The dashed line shows the expected /(oo) 
(scaled to 5xl0~^, roughly the expected abundance from 
stellar equilibrium chemistry, at very low mass-loss rates) 
for the parameters given in Sect.|Sl 

distance). The dependence on the luminosity is moderate. 
We therefore estimate that, within the adopted circum- 
stellar model, the derived SiO abundances are uncertain 
by at least a factor of three for those sources with multi- 
line observations. The uncertainty increases to a factor of 
five when only one transition is observed. 

For the high mass-loss rate (i.e., > 5x 10~® yr~^) 
Miras the situation is even worse. The radiation from these 
stars are significantly converted into longer-wavelength 
dust radiation, which has been taken care of only crudely 
by using two central blackbodies. Tests show that the re- 
sulting SiO line intensities are sensitive to the structure of 
the radiation sources, Sect 14.41 In addition, the SiO lines 
are rather saturated and hence the line intensities are, at 
least partly, insensitive to the abundance. Therefore, it is 
estimated that for these objects the SiO abundance is un- 
certain by a factor of five (in all cases information on three, 
or more, lines is available), but note that any reasonable 
change in the radiation field structure will systematically 
lower the abundance required to fit the data. 

7.3. Abundances 

It can be assumed that the stars in our samples have sili- 
con abundances clos e to th e solar value, Si/H = 3.6xl0~^ 
llAnders &: GrevesseL Il98l . If Si is fully associated with 
O as SiO, and all H is in H2, the maximum SiO fractional 
abundance is 7xl0~^. Detailed calculations on stellar at- 
mosphere equilibrium chemistry give abun dances in the 
vicinity of this for M-stars, about 4xl0~^ ijPuari et all 
Duari et al. also show that the SiO abundance is 
not affected by atmospheric shocks in the case of M-stars. 
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Table 4. Source parameters and SiO model results 



Source 


Var. 


P 


D 


[km s-"] 


«o(CO) 
[km s""] 


M 






/siO 


X?od 


N 




type 


[days] 


[pc] 


[10-^ Mq yr-"] 


[lO"'' cm] 


[lO"'' cm] 


[lo-''] 






RS And 


SRa 


136 


290" 


4.4 


4.4 


1.5 


2.2 


4.0 


16 




1 


UX And 


SRb 


400 


280" 


12.8 


12.8 


4 


5.9 


3.8 


12 




1 


e Aps 


SRb 


119 


110" 


4.0 


4.5 


0.4 


2.0 


2.2 


14 


0.3 


3 


TZ Aql 


Lb 




470" 


4.8 


4.8 


1 


2.2 


3.2 


15 


8.1 


2 


SV Aqr 


Lb 




470" 


8.0 


8.0 


3 


4.0 


4.5 


34 


2.1 


2 


T Ari 


SRa 


317 


310" 


2.4 


2.4 


0.4 


1.1 


2.8 


5.2 


2.9 


2 


RX Boo 


SRb 


340 


110" 


7.8 


9.3 


5 


4.8 


5.2 


8.0 


1.2 


6 


RV Cam 


SRb 


101 


350" 


5.8 


5.8 


2.5 


3.0 


4.6 


4.5 




1 


TX Cam 


M 


557 


380 


16.0 


18.5 


60 


33 


13 


5.5 


2.0 


3 


R Cas 


M 


431 


220 


7.0 


10.5 


13 


11 


7.4 


7.0 


0.8 


3 


UY Cet 


SRb 


440 


300" 


6.0 


6.0 


2.5 


3.1 


4.5 


6.0 


3.0 


3 


CW Cnc 


Lb 




280" 


7.0 


8.5 


5 


4.5 


5.1 


2.7 




1 


R Crt 


SRb 


160 


170" 


10.6 


10.6 


8 


6.1 


5.7 


6.0 


2.5 


3 


R Dor 


SRb 


338 


45" 


5.0 


6.0 


1.3 


3.3 


3.3 


5.0 


2.4 


4 


AH Dra 


SRb 


158 


340" 


6.4 


6.4 


0.8 


2.0 


2.6 


17 




1 


CS Dra 


Lb 




370" 


11.6 


11.6 


6 


5.0 


5.0 


2.7 




1 


S Dra 


SRb 


136 


270" 


9.6 


9.6 


4 


4.9 


4.4 


7.0 


4.3 


2 


SZ Dra 


Lb 




510" 


9.6 


9.6 


6 


5.0 


5.3 


1.8 




1 


TY Dra 


Lb 




430" 


9.0 


9.0 


6 


4.8 


5.6 


10 


0.3 


2 


R Hya 


M 


388 


150 


4.5 


7.0 


3 


4.0 


4.5 


7.0 


4.2 


2 


W Hya 


SRa 


361 


65" 


6.5 


6.5 


0.8 


2.8 


2.4 


15 


4.5 


7 


R Leo 


M 


313 


130 


6.0 


6.0 


2.0 


3.9 


5.5 


13 


2.5 


3 


U Men 


SRa 


407 


320" 


7.2 


7.2 


2.0 


3.4 


3.8 


5.8 


3.5 


2 


T Mic 


SRb 


347 


130" 


4.8 


4.8 


0.8 


2.2 


3.0 


5.3 


1.2 


2 


GX Mon 


M 


527 


540 


18.7 


18.7 


400 


48 


29 


0.8 


4.0 


5 


S Pav 


SRa 


381 


150" 


4.8 


9.0 


0.8 


3.8 


2.1 


2.6 


0.3 


2 


SV Peg 


SRb 


145 


190" 


6.3 


7.5 


3 


3.5 


6.5 


5.1 




1 


TW Peg 


SRb 


929 


200" 


9.5 


9.5 


2.5 


4.5 


3.6 


2.4 




1 


WX Psc 


M 


660 


600 


19.3 


19.3 


110 


22 


16 


6.0 


3.1 


4 


Pup 


SRb 


141 


85" 


2.3 


2.3 


0.2 


1.0 


2.1 


14 


2.1 


3 


Y Scl 


SRb 




330" 


5.2 


5.2 


1.3 


2.4 


3.6 


5.0 


0.7 


2 


V1943 Sgr 


Lb 




150" 


4.6 


5.4 


1.3 


2.8 


3.5 


7.3 


1.3 


2 


IK Tan 


M 


500 


250 


17.5 


18.5 


300 


31 


25 


0.4 


2.8 


4 


V Tel 


SRb 


125 


290" 


6.8 


6.8 


2.0 


3.2 


3.7 


5.0 


12.0 


2 


Y Tel 


Lb 




340" 


3.5 


3.5 


5 


1.6 


2.7 


54 


0.1 


2 


Y UMa 


SRb 


168 


220" 


4.8 


4.8 


1.5 


2.8 


4.1 


12 




1 


SU Vel 


SRb 


150 


250" 


5.5 


5.5 


2.0 


3.3 


4.3 


2.7 


5.3 


2 


BK Vir 


SRb 


150 


190" 


4.0 


4.0 


1.5 


1.8 


4.6 


2.3 




1 


RT Vir 


SRb 


155 


170" 


6.2 


7.8 


5 


3.6 


5.6 


13 


1.3 


2 


SW Vir 


SRb 


150 


120" 


7.5 


7.5 


4 


3.5 


4.8 


3.5 


9.1 


4 


IRC+10365 


M 


500 


750 


16.2 


16.2 


300 


46 


29 


4.0 


3.9 


3 


IRC- 10529 


M 


680 


270 


12.0 


12.0 


25 


6.8 


10 


1.1 


0.1 


3 


IRC-30398 


M 


575 


390 


16.0 


16.0 


60 


17 


13 


0.3 


10.0 


3 


IRC+40004 


M 


750 


410 


18.0 


18.0 


60 


19 


13 


0.2 


8.7 


3 


IRC+50137 


M 


629 


410 


17.0 


17.0 


100 


12 


15 


0.5 


2.0 


3 



" Distance derived assuming a luminosity of 4000 Lq 



The derived SiO abundances are given in Table 01 The 
distribution for the IRV/SRV sample has a median value 
of 6xl0~^, and a minimum of 2xl0~^ and a maximum 
of 5x10^^. For the IRVs and SRVs the median results are 
9x10^^ and 6x10^^, respectively. This is almost a factor 
of ten lower than expected from theory. Figure |31 shows 
the SiO abundance as a function of the mass-loss rate. In 
addition to the abundances being low, there is also a trend 



in the sense that both the upper and the lower 'envelope' 
of the abundances decrease with increasing mass-loss rate. 

The low mass-loss rate Miras follow the trend of 
the IRV/SRVs, and for the high mass-loss rate (M > 
5 X 10^^ M© yr~^) Miras we find a substantially lower 
abundance, a median below 10^^. Thus, the inclusion of 
the Miras shows that the trend of decreasing SiO abun- 
dance with increasing mass-loss rate continues towards 
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Fig. 6. Comparison a^tfee/^f^d CO (upper) and SiO ^(fS^fep/line profiles (in histograi]i'%[te/^r R Dor (left) and 
R Hya (middle) and GX Mon (right). The corresponding best- fit (i.e., to all observed line intensities) model SiO lines 
are also shown as solid lines 



high mass-loss rates. This is further discussed in Sect. |S1 
where an interpretation in terms of increased adsorption 
of SiO onto dust grains the higher the mass-loss rate is 
advocated. 

The spread in abundance, at a given mass-loss rate, is 
substantial, but it is within the (considerable) uncertain- 
ties, except possibly for the high mass-loss rate IVIiras, for 
which there seem to be a division into a low abundance 
group (on average 4x10"'') and a high abundance group 
(on average 5xl0~^), while . This division into two well- 
separated groups is peculiar, but within the circumstellar 
model used here this conclusion appears inescapable. One 
can argue that the modelling of the high mass-loss rate 
Mira SiO line emission is particularly difficult, but we find 
no reason why errors in the model should affect stars with 
essentially similar properties (L, M, Vc) so differently. 

7.4. CSE dynamics 

The observed SiO line profiles are used in the modelling 
to derive the gas expansion velocities in the regions of the 
CSEs where the observed SiO line emission stems from. A 
comparison of these values with the gas expansion veloc- 
ities derived from the modelling of circumstellar CO line 
emission is indeed a direct probe of the CSE dynamics 
since the extents of the SiO and CO line emitting regions 
are very different. 

The SiO and CO radio line profiles are clearly differ- 
ent, although this conclusion is mainly based on the lim- 
ited number of sources where the S /N-ratio of the data are 
high enough for both species. In Table ^ different values 



for the gas expansion velocity estimated from the SiO and 
the CO data are reported in the 11 cases where these are 
regarded as significantly different. In all cases the SiO ve- 
locities are smaller than those obtained from the CO data. 
Indeed, the SiO line profiles are narrower in the sense that 
the main fraction of the emission comes from a velocity 
range narrower than twice the expansion velocity deter- 
mined from the CO data. On the other hand, the SiO line 
profiles have weak wings so that the total velocity width 
of its emission is very similar to that of the CO emission. 
This is illustrated in Fig. El where we also show the cor- 
responding best-fit (i.e., to all observed line intensities) 
model SiO lines. It is clear that the model line profiles do 
not provide perfect fits to the observed line profiles, but 
they show that for the lower mass-loss rate objects the SiO 
line profiles are strongly affected by selfabsorption on the 
blue-shifted side. This explain partly why the SiO lines are 
narrower than the CO lines. The remaining discrepancy is 
interpreted as due to the influence of gas acceleration in 
the region which produces a significant f raction of the SiO 
l ine e mission, as suggested already by iBuiarrabal et alJ 
l)l986() . This interpretation is quantitatively corroborated 
by our modelling results when a velocity gradient is in- 
cluded, se e Sect. 14.41 The ext ent of the effect is though 
uncertain. teie ging et al.l l 200 dh. by comparing high- J SiO 
lines with CO line data, concluded that the SiO lines are 
formed predominantly in the part of the CSE where the 
gas velocity exceeds 90% of the terminal velocity. We sus- 
pect that the discrepancy between the widths of the SiO 
and CO lines decreases with the mass-loss rate of the ob- 
ject. In addition, we find that for at least some of the 
high-mass- loss-rate sources the higher- J SiO lines become 
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Fig. 7. The derived CO gas expansion velocities as a func- 
tion of the mass- loss rates for the IRV (squares), SRV (cir- 
cles), and Mira (triangles) samples. The dashed line shows 
the correlation found for the IRV/SRV sample (see text) 

essentially triangular, see GX Mon in Fig. O The model 
does a fairly good job in reproducing these SiO line pro- 
files, except that the model lines are less sharply peaked. 
A high sensitivity, multi-line study combined with inter- 
ferometric observations are required to fully tackle this 
problem. 

In this connection we also present Fig. [3 which shows 
the gas expansion velocity (determined from CO line mod- 
elling) as a function of mass loss rate for the IRV/SRV 
and Mira sa mples. This is an extension of the result of 
lOlofsson et al. 1.2002,) . and it shows that low to interme- 
diate mass- loss rate winds have a scaling of Vc oc 
and that this gradually goes over into a wind of close to 
20kms~^, for higher m ass-loss rates. This is a s expected 
for a dust-driven wind l)Elitzur &: Ivezi4 l200lj) . 

7.5. Peculiar sources 

To single out peculiar sources is a highly subjective pro- 
cess, and it also depends strongly on the S/N-ratio of the 
data (at high enough S/N-ratio probably most sources 
show a deviation from the expected). Here, a few sources 
in the IRV/SRV sample which qualify as peculiar or for 
which we have problems in the SiO line modelling are dis- 
cussed; 

iKerschbaum fc OlofssonI lll999l) found four objects in 
their sample of circumstellar CO radio line emission, which 
clearly show double-component line profiles, a narrow fea- 
ture centred on a broad plateau (E P Aqr, RV Boo, X Her , 
and SV Psc), aU of them SRVs. lOlofsson et all l)2002|) 
determined mass-loss rates and gas expansion velocities 
by simply decomposing the emission into two components 
and assuming that the emissions are additive. They found 
that the mass loss rates are higher for the broader com- 
ponent by, on average, an order of magnitude. The gas 



expansion velocities derived from the narrow components 
(wl.5 kms^^) put to question an interpretation in the 
form of a spherical o utflow. The origin of s uch a line profile 
is still not clear fsee lOlofsson et alTl(2002l) for a discussion 
on this issue). These four sources are also included in our 
SiO sample, and the spectra are shown in Figs. IB^ and 

ESI 

Towards EP Aqr there is no sign of the narrow fea- 
ture in the SiO J=2^1 and J=3^2 lines, only the broad 
feature is clearly present. This suggests that the broad fea- 
ture originates in a 'normal' CSE, while the narrow fea- 
ture may have a different origin. We note though that the 
SiO line profile of the broad component deviates somewhat 
from a smooth symmetric profile. SV Psc is very similar to 
EP Aqr in CO in the sense that the narrow feature is very 
much narrower than the broad feature. Unfortunately, the 
SV Psc SiO data are of low quality, but both components 
appear to be present. In the cases of RV Boo and X Her the 
CO and SiO line profiles are very similar, and the widths 
of the narrow components are about half of those of the 
broad ones. The SiO abundances of both components have 
been obtained, assuming that the emissions are additive. 
The results are given in Table [S] For all sources, and for 
both components, the results appear normal. 

Pup was singled out in lOlofsson et alJ l|2002l ) as a 
low mass-loss rate (2x 10~* Mq yr~^), low gas expansion 
velocity (2.1kms~ ^) object. This st ar has been recently 
discus sed also by lJura et aU l)2002j) and I Winters et alJ 
l|2002|) . In the latter paper comparisons are made with 
wind models, and it is concluded that stars with the mass- 
loss properties of Pup can be understood in terms 
of a pulsationally driven wind, where dust plays no dy- 
namic role. Our SiO line profiles resemble to some ex- 
tent those of CO in the sense that the narrow feature is 
also present. However, the SiO lines clearly show broad 
line wings. Fig. |H1 The full velocity width of these lines 
are wl2kms~^, i.e., larger than the CO line width, but 
narrower than the SiO(t'=l , J=2^1) maser lin e width 
of «20kms~"'^ measured bv lWinters et al.l l)2002(l . In ad- 
dition, the narrow feature, which appears narrower in 
the SiO lines than in the CO lines (Fig. IHJ, is not ex- 
actly centered on the broad component, its center lies at 
Vhc\ = 52.8kms~^ as opposed to 51.4kms~^ for the latter. 
This suggests a rather complicated dynamics in the inner 
part of the CSE, but high-quality data, also in higher-J 
SiO lines, are required before progress can be made. 

W Hya is one of the sources for which we have 
the highest quality data. It is also one of the sources 
with the poorest best-fit model. A much better fit is 
obtained by increasing the size of the SiO envelope to 
rc=6xl0^^ cm (and /sio=8xl0~^ as determined from the 
high- J lines), i.e., almost a factor of three higher than 
that obtained from Ea. l|13|l . Considering the uncertain- 
ties this is of no major concern. However, it is worth 
recalling that .Olofsson et al. (2002 ) derived a (molecu- 
lar hydrogen) mass- loss rate of 7x 10~^ Mq yr^^ from CO 
data [this result has been confirmed by including CO 
J=1^0 and 2^1 IRAM 30 m data ijBuiarrabal et all 
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Table 5. Source parameters and model results for those objects with double component line profiles 



Fig. S.L-^ Pup 

SiO 

line 

spec- 
tra 
and 
a 

CO(J=3^2) 
spec- 
trum 
(from 

lOlofsson et al.l 
ll2002D l 



Source 


Var. 

type 


P 

[days] 


D 
[pc] 


comp. 


M 

[10-^ M© yr-i] 


^;c(SiO) 
[kms-^] 


^;o(CO) 
[kms-i] 


[10^^ cm] 


./siO 

[io-«] 


N 


EP Aqr 


SRb 


55 


140 


broad 


5.0 


7.8 


9.2 


5.2 


3.8 


2 


RV Boo 


SRb 


137 


280 


broad 


2.0 


6.8 


7.0 


3.2 


6.0 


2 










narrow 


0.3 


3.0 


2.3 


2.6 


7.0 


1 


X Her 


SRb 


95 


140 


broad 


1.5 


6.5 


6.5 


3.4 


12 


1 










narrow 


0.4 


2.5 


2.2 


3.1 


4.0 


1 


SV Psc 


SRb 


102 


380 


broad 


3.0 


8.6 


9.5 


4.0 


10 


3 










narrow 


0.4 


1.6 


2.2 


3.2 


6.0 


1 



Il989t ICernicharo eTail ^9^, CO J=2^1, 3^2, and 
4^3 JCMT archiv e data, and the CO ISO result s of 
iBarlow et aO l(l99fi>) ]. while IZubko fc Elitzurl l(200(]l) re- 
quired a much higher mass-loss rate, 2.3x 10~^ Mq yr~^ 
(at the larger distance 115 pc) to explain the ISO H2O 
data. We have found that such a high mass-loss rate pro- 
duces CO radio lines that are at least a factor of 30 too 
strong. However, the ISO CO J=16^15 and J=17^16 
lines are only about a factor of two too strong. Hence, 
there is some considerable uncertainty in the properties 
of this CSE. A fit to the SiO line data using the larger 
distance and mass-loss rate is as bad as that for the low 
distance and mass-loss rate. 

In the case of R Dor lOlofsson et al.l l)2002j) could not 
fit well the CO radio line profiles. The model profiles 
were sharply double-peaked, while the observed ones were 
smoothly rounded. We merely note here that there was 



no problem to fit the SiO line profiles with the nominal 
values for R Dor. 



8. Discussion and conclusions 

An extensive radiative transfer analysis of circumstellar 
SiO 'thermal' radio line emission from a large sample of 
M-type AGB variable stars have been performed, partly 
based on a new, large, observational data base. It is con- 
cluded that, at this stage, the modelling of the circumstel- 
lar SiO radio line emission is considerably more uncertain 
than that of the CO radio line emission. Partly because 
the SiO line emission predominantly comes from the in- 
ner regions where the observational constraints are poor, 
but also partly because the behaviour of the SiO molecule 
is more complex, e.g., adsorption onto grains. A rather 
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detailed sensitivity analysis has been done, in order to es- 
timate the reliability of the derived results. 

In particular, the size of the SiO envelope is crucial to 
the modelling. Multi-line SiO modelling of eleven sources 
were used to establish a relation between the size of the 
SiO envelope and the density measure M jv^. This is of 
course rather uncertain, both in the absolute scale and 
in the dependence on the density measure. Comparison 
with estimates based on rather simple condensation and 
photodissociation theories suggests that the derived rela- 
tion is not unreasonable. A very good agreement with the 
photodissociation radii is obtained for an unshilded pho- 
todissociation rate of 2.5 x 10~^°s~^. It was also checked 
against interferometeric SiO line brightness size estimates 
of six sources. 

The SiO abundance distribution of the IRV/SRV sam- 
ple has a median value of 6xl0~^, and a minimum of 
2xl0~^ and a maximum of 5xl0~^. For these, low to 
intermediate mass-loss rate objects, we expect the abun- 
dances to be representative for the region inside the SiO 
photodissociation radius. This applies also to the low and 
intermediate mass-loss rate Miras. The high mass-loss rate 
Miras have a median abundance which is more than a fac- 
tor of six lower than that of the IRV/SRV sample. The 
derived SiO abundances are in all cases (within the uncer- 
tainties) below the abundance expected frorn stell ar at- 
mosphere chemistry f?^ 4x 10~^. iDuari et al.l l)l999|) l. the 
median for the total sample is down by a factor of ten. 
We regard this as a safe result, and interpret it in terms 
of SiO adsorption onto grains, which is efficient already at 
low mass-loss rates. 

In addition, there is a trend of decrasing SiO abun- 
dance with increasing mass-loss rate, Fig.|31 Here, we can- 
not entirely exclude systematic effects of the modelling. In 
particular, the adopted SiO envelope size relation can in- 
troduce such effects. E.g., smaller envelopes at low mass- 
loss rates, as indicated by the photodissociation model, 
would increase the estimated abundances for these ob- 
jects. This will actually strenghten the observed trend. 
There is no obvious reason for a similar envelope size de- 
crease for the high mass-loss rate objects (however, see 
below), but if present it would lead to a less pronounced 
trend. The discussion in Sect. 14.41 on the sensitivity on 
the radiation field distribution suggests that the abun- 
dance estimates for the high mass-loss rate objects are 
upper limits. Therefore, considering also that the effect 
is rather large, we regard the trend as at least tenta- 
tive. An interpretation in terms of increased adsorption 
of SiO onto grains with increasing mass-loss rate is nat- 
ural. In Fig. [S] a depletion curve based on the results in 
Sect. lfi.ll is plotted, and it represents well the general trend 
(the adopted parameters are Tbb = 2500K, L = 4000Lq, 
Uc = 10kms^^, 1/1 = 0.002, flgr = 0.05^m, pgr = 2gcm~^, 
Qp F = 0.03, rbi„d.siO = 29500 K, and a = l)t We empha- 
size once again how sensitive the theoretical condensation 
results are to the adopted parameters, and the depletion 
curve can easily be made to fit better the estimated abun- 
dances. 



For the high mass-loss rate Miras the SiO abundance 
distribution appears bimodal, a low abundance group (on 
average 4x10"'') and a high abundance group (on aver- 
age 5xl0~^). At this point we cannot identify any rea- 
son for this. The stars and their CSEs are rather similar 
and there is no reason to expect the modelling to arti- 
ficially produce very different results for rather similar 
objects, but the SiO line modelling of these objects are 
particularly difficult as discussed above. The high values 
can be explained if there is a process which decreases the 
condensation onto, or leads to effective evaporation from, 
dust grains for some objects. The former is possible if the 
dust-to-gas mass ratio is low. If, on the other hand, the 
dust-to-gas mass ratio is high, the region contributing to 
the SiO line emission may be much smaller than used in 
our modelling, and hence the SiO abundance is under- 
estimated. This may be the case for the low abundance 
objects. Substantial mass-loss rate variations with time 
may of course lead to surprising results. This can possi- 
bly be checked by high angular resolution observations of 
bot h CO and SiO radio line emission. It is also interesting 
that lWoods et al.l 1^2003^) found, in a sample of high mass- 
loss rate C-stars, that the SiO abundance is one of the 
few of their abundance estimates that vary significantly 
from star to star. We note, though, that for the C-stars 
the estim ated SiO abundances are lo w, about IxlO^*", 
and that IWillacv fc Cherchnefj \\9i9i ^ have shown that 
for C-stars shock chemistry may significantly alter the SiO 
abundance. The s a,me is not th e case for the 0-rich chem- 
istry according to 'Pua ri et al.l ljl999l) . but grains are not 
included in their analysis. 

The SiO and CO radio line profiles differ in shape. 
For those stars with high enough S/N-ratio data on both 
species, it is clear that the dominating parts of the SiO 
profiles are narrower than the CO profiles, but the for- 
mer have low-intensity wings which cover the full velocity 
range of the CO profile. The effect is more evident in high- 
J lines, and less evident in high mass-loss rate objects. 
This is interpreted (as has been done also by others) as 
due to the influence of gas acceleration in the region which 
produces most of the SiO line emission. This points to a 
weakness in our analysis. Clearly, this acceleration region 
must be treated more carefully in the radiative modelling, 
but this is also the region where condensation occurs, a 
process which is difficult to describe in detail. 

These results strongly suggest that SiO radio line emis- 
sion can be used as a sensitive probe of circumstellar dust 
formation and dynamics. However, considerable progress 
in this area can only be expected from a combination of 
high-quality SiO multi-line observations, high-quality in- 
terferometric observations of a number of SiO lines for 
a representative sample of sources, a detailed radiative 
transfer analysis, which includes also the dust radiation, 
and a detailed SiO chemical model. The rather high 
values of some of our best-fit models suggest that our cir- 
cui nstellar model needs t o be improved. 

lOlofsson et alJ l|2002f l identified a number of sources 
with peculiar CO line profiles, essentially consisting of 
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a narrow feature centered on a (much) broader feature. 
These have been discussed here from the point of view 
of their SiO line properties. Except in one case, the SiO 
and CO hne profiles are rather similar, and the derived 
SiO abundances are in no case peculiar. The low gas ex- 
pansion velocity source Pup has a very narrow SiO 
line profile as expected, but also a considerably broader, 
low-intensity component. Finally, W Hya imposes a prob- 
lem for the SiO line modelling. In principle, a much larger 
SiO envelope than warranted by the mass-loss rate derived 
from the CO data is required to fit well the (high-quality) 
SiO data. This, in combination with other data, suggest 
that the CSE of this star is not normal, possibly as an 
effect of time-variable mass loss. 
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Table A.l. Observational results of circumstellar SiO radio line emission towards a sample of M-type IRVs and SRVs 
(Part 1) 



GCVS4 


IRAS 


Code 


s 


/ 

[Kkms"^] 


-^mb 
L J 


^hcl 

[kms"^l 


'^LSR 

[kms"^] 


tic 

[kms"^] 


Q 


c 


BC And 


22586+4614 


021 


N 


0.2 










5 




RS And 


23528+4821 


021 


D 


0.71 


0.076 


-3.1 


4.2 


6.1 


3 




UX And 


02302+4525 


021 


D 


0.94 


0.044 


-19.4 


-20.1 


16.6 


3 




6 Aps 


14003-7633 


S21 


D 


0.62 


0.078 


10.8 


3.1 


6.1 


2 








S32 


D 


1.4 


0.18 


10.9 


3.2 


5.6 


2 








S54 


D 


2.7 


0.30 


10.7 


3.0 


7.6 


2 




TZ Aql 


20276-0455 


S21 


D 


0.15 


0.025 


49.3 


62.3 


5.0 


3 








S32 


D 


0.14 


0.026 


49.4 


62.4 


4.5 


4 




V584 Aql 


20079-0146 


S21 

S32 


N 

N 


0.2 

0.1 










5 
5 




AB Aqr 


22359 1417 


S21 
S32 


N 
N 


0.1 
0.1 










5 

5 




EP Aqr 


21439-0226 


S21 


D 


3.2 


0.28 


-41.3 


-31.9 


8.2 


2 








S32 


D 


4.4 


0.40 


-41.3 


-31.9 


7.9 


2 




SV Aqr 


23201-1105 


S21 


D 


0.54 


0.036 


7.3 


8.5 


9.5 


3 








S32 


D 


0.66 


0.070 


6.6 


7.8 


7.2 


2 




T Ari 


02455+1718 


021 


D 


0.084 


0.043 


6.5 


— 1.2 


1.5 


4 








S21 


D 


0.081 


0.021 


6.4 


-1.4 


2.9 


4 








S32 


D 


0.091 


0.017 


5.5 


-2.2 


4.0 


4 




RV Boo 


14371+3245 


021 


D 


0.79 


0.069 


-5.0 


9.5 


6.5 




b 






021 


D 


0.19 


0.051 


—6.0 


8.5 


2.0 




n 






021 


D 


0.98 


0.12 


-6.9 


7.6 


8.8 


2 


b+n 






132 


D 


2.5 


0.25 


-5.9 


8.6 


8.0 




b+n 


RX Boo 


14219+2555 


021 


D 


9.5 


0.84 


-12.1 


0.2 


9.0 


1 








S21 


D 


5.1 


0.44 


-11.6 


0.7 


8.7 


1 








S32 


D 


6.4 


0.61 


-11.7 


0.6 


8.2 


1 




RV Cam 


04265+5718 


021 


D 


0.32 


0.036 


-15.9 


-16.0 


6.8 


4 




BI Car 


10416-6313 


S21 
S32 


N 
N 


0.2 
0.3 










5 
5 




SS Cep 


03415+8010 


021 


N 


0.3 










5 




V744 Cen 


13368-4941 


S21 
S32 


N 
N 


0.1 
0.4 










5 
5 




V806 Cen 


13465-3412 


S21 
S32 


N 
N 


0.3 
0.2 










5 

5 




UY Cet 


00245-0652 


021 


D 


0.41 


0.071 


6.5 


3.6 


4.6 


4 








S21 


D 


0.38 


0.053 


7.8 


4.9 


5.2 


4 








S32 


D 


1.02 


0.11 


8.9 


6.0 


7.1 


3 




CW Cnc 


09057+1325 


021 


D 


0.62 


0.065 


26.6 


17.2 


7.9 


3 




R Crt 


10580-1803 


S21 


D 


5.0 


0.33 


19.6 


11.9 


10.9 


1 








S54 


D 


9.9 


0.66 


19.8 


12.1 


10.7 


2 








S65 


D 


13.6 


0.86 


19.0 


11.3 


12.6 


3 




V CVn 


13172+4547 


021 


N 


0.3 










5 




W CvK 


21341+4508 


021 


N 


0.1 










5 




U Del 


20431+1754 


021 


N 


0.5 










5 




R Dor 


04361-6210 


S21 


D 


9.6 


1.1 


23.9 


7.5 


6.1 


1 








S32 


D 


24.0 


2.8 


24.0 


7.6 


6.0 


1 








S54 


D 


31.1 


3.6 


24.4 


8.0 


6.0 


1 








S65 


D 


34.7 


5.1 


24.4 


8.0 


5.8 


1 




AH Dra 


16473+5753 


132 


D 


1.2 


0.17 


57.6 


74.9 


8.5 






CS Dra 


11125+7524 


021 


D 


0.41 


0.022 


-60.1 


-51.3 


12.5 


4 




S Dra 


16418+5459 


021 


D 


1.4 


0.10 


-2.2 


15.4 


9.9 


4 








132 


D 


3.6 


0.30 


— 1.6 


15.9 


8.6 


3 




SZ Dra 


19089+6601 


021 


D 


0.17 


0.024 


-44.2 


-28.3 


4.8 


4 




TY Dra 


17361+5746 


021 


D 


0.65 


0.038 


-33.9 


-18.9 


11.5 


4 








132 


D 


3.2 


0.22 


-34.0 


-16.4 


11.1 


3 




UU Dra 


20248+7505 


021 


N 


0.3 










5 




g Her 


16269+4159 


021 


N 


0.2 










5 




X Her 


16011+4722 


021 


D 


2.2 


0.18 


-90.5 


-73.0 


8.0 




b 






021 


D 


0.53 


0.16 


-89.5 


-72.0 


2.2 




n 






021 


D 


2.7 


0.34 


-90.7 


-73.2 


7.4 


1 


b+n 


FZ Hya 


08189+0507 


021 


N 


0.7 










5 
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Fig. B.l. SiO spectra of M-type IRVs. Note tiie heliocentric velocity scale 



20 D. Gonzalez Delgado et al.: SiO thermal emission from M-type AGB stars 

Table A. 2. Observational results of circumstellar SiO radio line emission towards a sample of M-type IRVs and SRVs 
(Part 2) 



GCVS4 


IRAS 


Code 


s 


7 

[Kkms"^] 


L J 


ViiqI 

[kms"^] 


[kms"^] 


Vc 

[kms"^] 


Q 


W Hya 


13462-2807 


S21 


D 


6.0 


0.63 


38.5 


40.4 


7.1 


1 






S32 


D 


10.9 


1.2 


39.2 


41.1 


6.4 


1 






S54 


D 


13.2 


0.98 


39.2 


41.1 


7.1 


1 






S65 


D 


13.2 


0.86 


39.4 


41.3 


6.7 


1 


RW Lep 


05365-1404 


S21 


N 


0.3 










5 






S32 


N 


0.1 










5 


U Men 


04140-8158 


821 


D 


0.30 


0.037 


28.9 


17.0 


6.3 


3 






S32 


D 


0.44 


0.049 


28.7 


16.8 


6.5 


3 


T Mic 


20248-2825 


S21 


D 


0.67 


0.088 


17.7 


25.3 


5.6 


2 






S32 


D 


1.2 


0.16 


17.8 


25.4 


5.6 


2 


EX Ori 


05220 0611 


021 


N 


0.7 










5 


V352 Ori 


05592-0221 


021 


N 


0.5 










5 


S Pav 


19510-5919 


S21 


D 


0.32 


0.030 


-19.2 


-20.0 


8.3 


3 






S32 


D 


0.40 


0.055 


-19.1 


-19.9 


5.6 


4 


NU Pav 


19575-5930 


S21 


N 


0.1 










5 






S32 


N 


0.3 










5 


SV Peg 


22035+3506 


021 


D 


2.1 


0.18 


-8.3 


2.1 


9.5 


2 


TW Pee 


22017+2806 


021 


D 


0.54 


0.035 


—23.0 


—10.5 


9.3 


4 


SV Psc 


01438+1850 


021 


D 


0.45 


0.05 


9.5 


5.7 


6.9 


3 






S21 


D 


0.31 


0.03 


10.5 


6.7 


11.0 


3 






S54 


D 


0.63 


0.06 


10.4 


6.6 


8.4 


3 






S54 


D 


0.20 


0.10 


10.4 


6.6 


1.6 


3 






S54 


D 


0.83 


0.10 


10.2 


6.4 


6.3 


3 


V PsA 


22525-2952 


S21 


N 


0.1 










5 






S32 


N 


0.1 










5 


Pup 


07120-4433 


S21 


D 


0.83 


0.14 


52.4 


33.6 


4.6 


2 






S32 


D 


1.7 


0.30 


52.0 


33.3 


3.5 


2 






S54 


D 


2.5 


0.54 


52.3 


33.6 


3.6 


1 


Y Scl 


23063-3024 


S21 


D 


0.22 


0.024 


30.9 


29.5 


6.4 


4 






S32 


D 


0.28 


0.036 


31.8 


30.4 


5.8 


3 


CZ Ser 


18347-0241 


021 


N 


0.4 










5 


r^Ser 


15341+1515 


021 


N 


0.2 










5 






S21 


N 


0.2 










5 






S32 


N 


0.1 










5 


SU Sgr 


19007-2247 


S21 


N 


0.1 










5 






S32 


N 


0.2 










5 


V1943 Sgr 


20038-2722 


S21 


D 


0.96 


0.13 


-23.1 


-14.4 


5.4 


2 






S32 


D 


1.7 


0.26 


-23.0 


-14.3 


4.8 


2 


V Tel 


19143-5032 


S21 


D 


0.59 


0.063 


-34.6 


-31.8 


7.1 


3 






S32 


D 


0.40 


0.060 


-35.0 


-32.2 


5.4 


3 


Y Tel 


20165-5051 


S21 


D 


0.21 


0.034 


—46.7 


—45.3 


3.8 


4 






S32 


D 


0.44 


0.066 


-45.8 


-44.4 


4.6 


3 


A 7 TIMa 




021 


N 


0.2 










5 


Y UMa 


12380+5607 


021 


D 


1.2 


0.14 


8.8 


18.6 


6.0 


2 


SU Vel 


09480-4147 


821 


D 


0.35 


0.035 


20.8 


7.0 


7.4 


3 






832 


D 


0.38 


0.041 


20.8 


7.0 


7.4 


3 


BK Vir 


12277+0441 


021 


D 


0.63 


0.10 


15.1 


17.8 


4.7 


3 


RT Vir 


13001+0527 


021 


D 


3.8 


0.39 


13.5 


18.7 


7.4 


2 






821 


D 


3.2 


0.27 


13.4 


18.6 


9.2 


2 


RW Vir 


12046-0629 


021 


N 


0.4 










5 


SW Vir 


13114-0232 


021 


D 


3.5 


0.33 


-14.5 


-9.8 


7.9 


2 
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Table A. 3. Observational results of circumstellar CO and SiO line emission towards our Mira sample 
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GCVS4 


IRAS 


Molecule 


Gode 


I 
















[K kms"^] 


[K] 


[kms^^] 




TX Cam 


04566+5606 


GO 


OlO 


21.5 


0.88 


11.5 


20.0 








no 


67.0 


1.44 


11.5 


21.0 








121 


241.3 


3.93 


11.4 


20.7 








J21 


69.7 


2.61 


11.4 


20.2 








J32 


167.2 


6.14 


11.3 


20.3 






SiO 


021 


13.2 


0.52 


10.7 


18.6 


R Gas 


23558+5106 


GO 


OlO 


8.4 


0.46 


24.9 


12.2 








J21 


32.1 


1.78 


24.8 


11.1 








J32 


100.2 


5.48 


25.2 


11.0 








J43 


89.6 


5.44 


24.8 


11.9 






SiO 


021 


8.6 


0.67 


26.3 


9.4 


R Hya 


13269-2301 


GO 


SIO 


0.6 


0.09: 


-9.7: 


5.3: 








J32 


43.0 


3.65 


-10.5 


8.6 






SiO 


S21 


2.8 


0.44 


-11.5 


5.2 








S54 


7.7 


1.30 


-11.2 


5.9 


R Leo 


09448+1139 


GO 


OlO 


2.4 


0.25 


0.2 


7.4 








J21 


15.0 


1.30 


-0.4 


8.7 








J32 


41.6 


3.88 


-0.3 


7.8 






SiO 


021 


5.1 


0.51 


0.5 


6.8 








S21 


4.9 


0.65 


0.2 


6.9 








S54 


10.8 


1.58 


0.4 


6.3 


GX Mon 


06500+0829 


GO 


OlO 


49.6 


1.55 


-9.3 


24.4 








J21 


67.2 


2.27 


-9.5 


18.2 








J32 


87.6 


2.96 


-9.4 


17.9 








J43 


56.2 


1.82 


-8.7 


17.5 






SiO 


021 


9.4 


0.36 


-9.6 


19.5 








S21 


5.5 


0.29 


-9.6 


18.9 








S54 


8.9 


0.42 


-9.1 


18.7 


WX Psc 


01037+1219 


GO 


OlO 


52.0 


1.61 


9.6 


18.6 








J21 


31.0 


1.31 


10.1 


18.3 








J32 


45.5 


1.63 


9.5 


20.8 








J43 


49.9 


1.82 


9.5 


20.6 






SiO 


021 


10.0 


0.40 


10.1 


18.9 








S21 


5.6 


0.22 


9.3 


19.2 








S54 


6.8 


0.29 


9.3 


17.7 


IK Tau 


03507+1115 


GO 


OlO 


58.6 


1.70 


34.4 


17.2 








J21 


103.8 


3.32 


34.5 


17.5 








J32 


143.8 


4.81 


34.2 


17.6 








J43 


127.0 


4.52 


33.2 


18.0 






SiO 


021 


16.1 


0.66 


33.9 


18.1 








S54 


18.9 


0.73 


34.7 


17.5 








S65 


13.7 


0.61 


34.7 


15.9 


IRC+ 10365 


18349+1023 


GO 


OlO 


21.5 


0.79 


-31.3 


16.9 








,12 1 


42.0 


1.48 


-30.7 


15.3 






SiO 


021 


6.3 


0.30 


-32.6 


15.6 


IRC-10529 


20077-0625 


GO 


OlO 


15.7 


0.65 


-18.1 


15.7 








J21 


24.4 


1.18 


-17.5 


17.7 








J32 


48.0 


2.22 


-17.1 


14.2 








J43 


57.4 


2.87 


-17.2 


16.0 






SiO 


021 


2.9 


0.16 


-17.7 


14.6 


IRC-30398 


18560-2954 


GO 


J21 


44.6 


1.76 


-6.4 


19.4 






SiO 


S21 


0.8 


0.06: 


-10.1: 


15.3: 


IRC+40004 


00042+4248 


GO 


OlO 


24.8 


0.82 


-20.5 


19.2 








J21 


41.5 


1.41 


-20.5 


18.9 






SiO 


021 


1.3 


0.05 


-21.3 


20.3 


IRC+50137 


05073+5248 


GO 


J21 


36.8 


1.40 


3.5 


18.7 








J32 


35.7 


1.32 


3.3 


18.1 






SiO 


021 


2.0 


0.10 


1.4 


14.4 
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Fig. B.2. SiO spectra of M-type SRVs (Part 1). Note tlie heliocentric velocity scale 
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Fig. B.3. SiO spectra of M-type SRVs (Part 2). Note tlie heliocentric velocity scale 
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Fig. B.4. SiO spectra of M-type SRVs (Part 3). Note tlie heliocentric velocity scale 
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Fig. B.5. SiO spectra of M-type Miras. Note the LSR velocity scale 



